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First-principles calculations of the electronic structure and vibrational modes, in 
a system of graphene bilayers and trilayers intercalated with alkaline earth atoms, 



O are presented. It is found that, in similarity to the case of superconducting graphite 

^ intercalation compounds, the Fermi level is crossed by an s-band derived from the 

^ intercalant states, as well as graphitic 7r-bands. The electron-phonon coupling pa- 

^ rameter A is found to be 0.60 and 0.80, respectively, in calcium intercalated graphene 

5^ bilayers and trilayers. In superconducting CaCg graphite intercalation compound, 

H the calculated value for A is 0.83. It is concluded that two-dimensional superconduc- 

_JL tivity is possible in a system of a few graphene layers intercalated with calcium. 
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Keywords: graphene intercalates, electron-phonon coupling, superconductivity 



> 

in 

^ I. INTRODUCTION 

1^ Superconductivity in graphite intercalation compounds (GICs) has been stud- 

O ied extensively for a long period of time [1-10]. A fascinating aspect of these 

, superconducting compounds is that they are usually formed from components 

;;• neither one of which is a superconductor. In the case of alkali metal GICs the 

• ^ superconducting transition temperature is of the order of 1 K. Within the GIG 

r\ Brillouin zone (BZ) there exists an alkali-derived spherical band centered at the F- 

d point and carbon-derived cylindrical vr-bands centered along the HK axis. It was 

proposed [6, 11, 12] that the presence of both bands is necessary for superconduc- 
tivity, and that a two-band model is capable of explaining various experimental 
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observations relating to the occurrence of superconductivity in alkali-metal GICs 
and their anisotropic magnetic properties. 

It was also found that alkaline earth GICs exhibit superconductivity [7, 8] with 
a transition temperature of 6.5 K (YbCe) and 11.5 K (CaCe). Careful analysis 
[13-16] of the electron-phonon coupling in CaCe reveals that this superconductor 
is of the conventional BCS type, and that the presence of both Ca-derived and 
C-derived bands is necessary to explain the relatively high superconducting tran- 
sition temperature in this compound. It should be pointed out, however, that 
some controversy still exists regarding the mechanism for superconductivity in 
CaCg GIC, whereas a large Ca isotope effect was measured in CaCg GIC [17, 18] 
indicating a dominant role played by phonons due to Ca atomic vibrations, angle- 
resolved photoemission spectroscopy (ARPES) measurements on CaCe GIC [19], 
on the other hand, indicate that graphitic high frequency phonon modes are most 
strongly coupled to electrons. 

The recent discovery of stable graphene sheets [20] raises the question of 
whether two-dimensional superconductivity is possible in intercalated graphene 
layers. The possibility of intercalating a system consisting of a few graphene 
layers with intercalants, such as FeCla and Br, has been recently demonstrated 
[21, 22]. It is clear that a system consisting of many graphene layers, if inter- 
calated, should exhibit properties similar to those of GICs. It is natural to ask 
whether a system consisting of a few graphene layers, intercalated with alkali, 
or alkaline earth, metal atoms, will exhibit superconductivity. This problem was 
discussed recently by Mazin and Balatsky [23] , who considered a graphene bilayer 
intercalated with Ca. They argued, on the basis of similarities with the electronic 
energy bands of the CaCe GIC, that a graphene bilayer intercalated with Ca is 
probably a superconductor. 

In this paper we examine the problem of superconductivity in intercalated 
graphene layers. In particular, we calculate the energy bands in a graphene bi- 
layer intercalated with the alkaline earth atoms Ca, Sr, and Ba. In all three cases 
we find that a band, derived from the alkaline earth atoms, crosses the Fermi en- 
ergy, along with carbon-derived vr-bands. We then calculate the electron-phonon 
coupling strength A in these systems. The values obtained are reasonably large; 
in the case of Ca-intercalated bilayer, we find A = 0.60, a value to be compared 
with 0.83 calculated in CaCg GIC [24]. Next, the case of a graphene trilayer 
intercalated with Ca is examined; here we find a value of 0.80 for the electron- 
phonon coupling strength A. These calculations indicate that superconductivity 
is indeed possible in intercalated graphene layers. 

In Section II we discuss the computational methods employed in this work. 
The results of the energy bands, vibrational modes, and electron-phonon coupling 
calculations are presented and discussed in Section III, and conclusions are given 
in Section IV. 



II. METHODS 

The total energy and electronic band structure calculations are carried out 
using the all-electron, full-potential, linear augmented plane wave (FP-LAPW) 
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method as implemented in the WIEN2K code [25]. The exchange-correlation 
potential was calculated using the generalized gradient approximation (GGA) as 
proposed by Pedrew, Burke, and Ernzerhof (PBE) [26]. The radii of the muffin-tin 
spheres for the carbon and the alkaline earth atoms were taken as 1.3ao and 2. Sag, 
respectively, where oq is the Bohr radius. We set the parameter RurKmax = 7, 
where Rmt is the smallest muffin-tin radius and Kmax is a cutoff wave vector. 
The valence electrons wave functions inside the muffin-tin spheres are expanded 
in terms of spherical harmonics up to Imax = 10, and in terms of plane waves 
with a wave vector cutoff K^ax in the interstitial region. The charge density is 
Fourier expanded up to a maximum wave vector Gmax = ISa^^. Convergence of 
the self-consistent field calculations is attained with a total energy convergence 
tolerance of 0.01 mRy. 

The calculations of the frequencies of the vibrational modes and the electron- 
phonon coupling parameter are performed using density functional theory (DFT) 
and ultrasoft pseudopotentials [27, 28]. The electron-phonon coupling strengths 
are calculated using density functional perturbation theory within the linear re- 
sponse approximation. Ultrasoft pseudopotentials [29] are used for both, the 
alkaline earth and the carbon atoms. The valence electrons wave function and 
charge density are expanded in plane waves using 30 Ry and 300 Ry cutoffs, 
respectively. In the calculation of the vibrational frequencies the electronic in- 
tegration is carried out using a uniform mesh of Nk = 8x8 k-points in the 
two-dimensional BZ. On the other hand, the electronic density of states, used in 
computing the electron-phonon coupling parameter A, is calculated using a finer 
uniform mesh of 24 x 24 k-points. The value of A is obtained by averaging over 
a uniform mesh of A^g = 6 x 6 phonon momentum q-points. 

III. RESULTS AND DISCUSSION 

The ffist issue to decide is this: given a graphene bilayer, will it be energetically 
favorable for the alkaline earth atoms to be intercalated between the two carbon 
sheets comprising the bilayer? Or will the energy be lower if the alkaline earth 
atoms sit on top of the bilayer? If the atoms sit on top of the bilayer, the resulting 
structure will be denoted by Ce-Cg-A, where A represents an alkaline earth atom. 
On the other hand, the intercalated structure is denoted by Cg-A-Cg. In order 
to decide which structure is more stable, we carried out total energy calculations 
for both structures. In the case of Cg-Cg-A, the alkaline earth atom is placed 
over the center of a hexagon in the upper sheet, and its distance to the sheet is 
varied in order to find the lowest energy. On other hand, in the case of Cg-A- 
Cg, the distance between the carbon sheets is fixed at the measured separation 
in the corresponding graphite intercalation compound [9, 30]. In all the cases 
considered, corresponding to A=Ca, Sr, and Ba, we found that the energy is lower 
for the intercalated structure. For the case A=Ca, the intercalated structure Cg- 
Ca-Cg is lower in energy than the Cg-Cg-Ca structure by 0.47 eV/Ca atom. In the 
other two cases, A=Sr and Ba, the energy of the Cg-A-Cg structure is lower than 
the Cg-Cg-A structure by 0.43 e\^/Sr atom and 0.59 e\^/Ba atom. We conclude 
that it is energetically favorable for the alkaline earth atom to be intercalated in 



between the graphene layers. 

(a) M^^^ C"' ^ 






1 


/ 


v 


1 


; < 


f 

• 


> 


(iM' 


1 > 

1 


\ ' 


/ 


1 


1. 




J 




FIG. 1: (Color online) (a) The unit cell of the structure Cg-A-Ce, where A is an alkaline 
earth atom, is shown, (b) The 2-D Brillouin zone of a graphene layer (or bilayer) is 
shown as a dashed hexagon, while that for the intercalated graphene bilayer is shown 
as a solid hexagon. Note that the M and K points in the graphene BZ (here denoted 
by M' and K') are folded onto points M and F, respectively, of the intercalated bilayer 
BZ. 

In the following we focus attention only on the Ce-A-Ce structure. The sepa- 
ration between the two carbon sheets of the intercalated bilayer is taken as 4.52 
A, 4.94 A, and 5.25 A for A=Ca, Sr, and Ba, respectively. The unit cell is a 
(•\/3a X \/3a) R 30° supercell of the graphene sheet unit cell; i.e., the in-plane 
lattice constant of Cq-A-Cq is a/3 times the lattice constant in the graphene sheet, 
and the in-plane lattice vectors of Ce-A-Ce are rotated by 30° relative to those 
in the graphene sheet. The intercalated bilayer unit cell, along with the first 
Brillouin zones of the graphene sheet and the supercell structure are shown in 
Figure 1. Since the DFT calculations are carried out on a 3-diniensional crystal, 
the c-lattice constant is taken to be sufficiently large (c = 14.6 A, 15.0 A, and 
15.3 A for A= Ca, Sr, and Ba, respectively) so that adjacent intercalated bilayers 
have negligible interactions. 

In Figure 2 the calculated energy bands in Cg-Ca-Cg are plotted in the fatband 
representation in two different panels, along high symmetry directions in the first 
BZ. The energy band dispersion curves in both panels are identical, but the width 
(or the "fatness") of the band at every k-point in Figure 2a is proportional to 
the contribution of the Ca atoms to the band eigenstates, while in Figure 2b, the 
width is proportional to the contribution of the C-derived states. We note that 
an s-band originating from the Ca states, as well as vr-bands originating from C 
states, cross the Fermi level. This is indicative of the partial charge transfer from 
the intercalant atoms to the carbon sheets. Indeed, using the Bader's atoms in 
molecule (AIM) method[31], we calculated a charge transfer to the carbon sheets 
of 1.2 electrons/Ca atom, 1.06 electrons/Sr atom, and 0.91 electrons/Ba atom. 
We also note that in the (v3 x v3) R 30° supercell structure, the K-point of the 
BZ of the graphene sheet is folded onto the F point of the BZ of the supercell, 
as Figure 1 indicates. This explains the observation that the vr-bands that are 
centered at the K-point in graphene are now centered at the F-point in the band 
structure shown in Figure 2. The total and atomic densities of states are plotted 
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FIG. 2: (Color online) Plot of the energy bands in Ca-intercalated graphene bilayer 
showing the character of the contributing states. The width of the band at every k- 
point is proportional to the contribution of (a) Ca-derived states, and (b) C-derived 
states. 



in Figure 3. We note that the total density of states (DOS) consists of the sum 
of the atomic DOS and the DOS in the interstitial region. The plot in Figure 3 
reveals that there is a substantial contribution to the total DOS coming from the 
interstitial region. 

It has been argued previously that for superconductivity in GICs to exist, 
it is necessary that the charge transfer from the intercalant to carbon sheets 
be incomplete [11, 12], so that a band derived from the intercalant states must 
coexist with the vr-bands derived from the carbon states. Superconductivity was 
then explained as resulting from scattering of electrons, through phonon emission 
or absorption, between the intercalant and the vr-bands. More recently, DFT 
calculations showed that in the CaCe GIG, the vr-bands alone can not account for 
the observed superconducting transition temperature, and that the presence of 
the Ga-derived band is necessary [13-16]. Hence, it is reasonable to expect that 
superconductivity is possible in the intercalated bilayers under discussion in this 
work. The superconducting transition temperature will be determined mainly by 
the strength of the electron-phonon coupling in these systems. 

Similar to the case of GaGg GIGs [24], we calculate the phonon spectrum of the 
bilayers using DFT as implemented in the PWSGF code [27, 28]. The technical 
details of the calculation are given in Section II. 

For the graphene bilayer intercalated with alkaline earth atoms, the crystal 
point group is D^h = Dq x i, where i is the inversion symmetry. The group's 
character table is given in Table I. The Cq axis is the 2;-axis which passes through 
the alkaline earth atom and is perpendicular to the carbon sheets. The Cg axes 
in the x — y plane are parallel to the lines bisecting the hexagonal sides, and the 
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FIG. 3: (Color online) Total and atomic densities of states in Ca-intercalated graphene 
bilayer. The total density of states is the sum of the atomic densities of states and the 
density of states in the interstitial region. The Fermi energy is at E= 0. 



C2 axes are parallel to the hexagonal diagonals. The 3A^ = 39-diniensional rep- 
resentation, r, is given at the bottom of Table I; it is obtained by considering the 
transformation properties of the 39 cartesian unit vectors placed on the thirteen 
atoms of the unit cell. From Table I, the 39 phonon modes at the BZ center are 
decomposed into irreducible representations according to 



phonon 



2Au + A2a + 2Bu + B2a + 3Ei, + ?,E. 



129 



'19 



29 



^19 



^29 



+ Alu + Blu + ?>A2u + 2B2U + 4Ei„ + 3E; 



2u 



The Aig, Eig, and E2g modes are Raman-active while the A2U and Eiu modes 
are infrared-active. From the representations of x, y, and z, in Table I, the 
three translational modes are found to transform as A2U and Ei^. The Raman 
and infrared-active modes, with nonzero frequencies, for the three intercalated 
bilayers considered in this work are given in Table II. The overall trend displayed 
in Table II is in accordance with what is expected on general grounds. The high 
frequency modes are lower than the corresponding modes in graphite as a result 
of the considerable charge transfer which softens the force constants in the carbon 
sheets. The frequencies of the higher modes increase slightly as the intercalant 
is changed from Ca to Sr to Ba, a consequence of the concomitant decrease in 
the amount of charge transfer. The validity of the first-principles calculations 
of the phonon spectra in the intercalated bilayer may be checked by repeating 
the same calculation for pristine graphite where the F point phonon frequencies 
are known experimentally. In graphite, the calculated phonon frequencies of the 
Raman-active modes are oj{E2g) = 75 cm~^ and 1557 cm~^, the IR-active modes 
<^{A2u) = 892 cm~^, u{Eiu) = 1567 cm~^, and the silent modes u{Big) = 133 



TABLE I: Character table of D^h = Dq x i. The character of the 39-dimensional 
representation, T, is given at the bottom. 



Axes 


Modes 


B 


C2 


2C3 


2Ce 


3^2 


3^2 


i 


iCs 


2iC3 


2iCe 


3iC^ 


siOi 




Alg 




1 
























Aiu 




1 










-1 


-1 


-1 


-1 


-1 


-1 




Mg 




1 






-1 


-1 










-1 


-1 


z 


Mu 




1 






-1 


-1 


-1 


-1 


-1 


-1 








Big 




-1 




-1 




-1 




-1 




-1 




-1 




Blu 




-1 




-1 




-1 


-1 




-1 




-1 






B2g 




-1 




-1 


-1 






-1 




-1 


-1 






B2u 




-1 




-1 


-1 




-1 




-1 






-1 




Big 


2 


-2 


-1 
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-2 


-1 










{x,y) 


Blu 
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-2 


-1 










-2 


2 




-1 










B2g 
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-1 


-1 
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-1 


-1 










B2u 
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-1 


-1 








-2 


-2 
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39 


-1 
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-1 


-1 


-3 


1 





-2 
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1 



crn"^ and 889 cm^^. Experimentally, the measured frequencies [32-35] of the 
Raman-active modes are 42 and 1582 cm~^, the IR-active modes 868 and 1588 
cm~^, and the low frequency silent mode 128 cm~^. It is seen that, except for 
the low frequency mode at 42 cm~^, a good agreement exists between the first- 
principles calculations and experiment, the error being of the order of a few 
percent. 

TABLE II: Raman- and Infrared-active modes in the intercalated graphene bilayer 
Ce-A-Ce for A=Ca, Sr, and Ba. 



a;(irrep)/cm 



-1 



Raman-active modes 


Ca 


69(^lg) 


l28{Alg) 


497(^1^) 


498(^2,;) 






1132(^23) 


1135(^lg) 


1343(^13) 


1410(£;23) 




Sr 


MBlg) 


133(Aig) 


506{Eig) 


511(^25) 






1152(^25) 


n52{Eig) 


1355(^13) 


1413(^23) 




Ba 


^8{Blg) 


115(Aig) 


506(^13) 


509(^2^) 






1166(^lg) 


1167(^29) 


1369(^13) 


1454(^23) 


Infrared-active modes 


Ca 


115(Si„) 
1349(A2„) 


293(^2«) 


481(^1,) 


1130(^i„) 




Sr 


72(^i„) 
1358(^2^) 


207(^2«) 


498 (Ei„) 


1151(^i„) 




Ba 


G7iBiu) 
1370(A2„) 


153(^2j 


503(^1^) 


1165(^i„) 



We calculated the electron-phonon coupling parameter A in Ce-A-Ce and found 
it to be 0.60, 0.42, and 0.35, for A=Ca, Sr, and Ba respectively. The transition 



superconducting temperature Tc of the the studied compounds was estimated 
using the Allen-Dynes [36] formula with values for the phonon logarithmic average 
frequencies uiog = 29.4 meV, 30.1 meV and 28.4 meV for Ce-Ca-Ce, Ce-Sr-Ce, 
and Ce-Ba-Ce respectively. Using a screened Coulomb pseudopotential /i* = 0.14 
we obtained Tc = 4.7 K for Ca-intercalated graphene bilayer, 0.95 K for Sr- 
intercalated bilayer, and 0.1 K for Ba- intercalated bilayer. This is to be compared 
with the value A = 0.83 calculated previously in CaCe GIC [24] for which Tc is 
11.5 K. 

In order to study the effect of increasing number of graphene layers on the 
superconducting properties, we consider a trilayer of graphene intercalated with 
calcium: Ca2Ci8. The unit cell is again given by {\/3a x y/3a) R 30° supercell of 
the graphene sheet unit cell. In conformity with the corresponding situation in 
stage-1 GICs, the carbon layer stacking is taken to be AAA. Between the first and 
second carbon layers, Ca atoms sit above hexagon centers, occupying positions 
labeled a. Between the second and third carbon layers, Ca atoms sit above 
hexagon centers, occupying positions labeled (3. The two intercalant layers are 
not directly above each other; the second intercalant layer is shifted, relative to 
the first intercalant layer, by a lattice vector of the graphene sheet, as illustrated 
in Figure 4. Thus, the stacking in the intercalated graphene trilayer is AaA(]A. 
The crystal point group is reduced, from Dq^ in the case of intercalated bilayer, 
to D^d in the case of the intercalated trilayer. The phonon modes, at the BZ 
center, are decomposed into the following irreducible representations 

Tphonon = 5Aig + 5A2g + AA^^ + QA2U + 10^3 + 10E„. 




FIG. 4: (Color online) Projection, on the plane of the first carbon layer, of the in- 
tercalated graphene trilayer. In the first intercalant layer, Ca atoms occupy positions 
denoted by a, whereas Ca atoms occupy positions /3 in the second intercalant layer. 

The Aig and Eg modes are Raman-active, whereas A2u and Eu modes are 
IR-active. The three zero-frequency translational modes transform according to 
A2u and E^ irreducible representations. In Table III we report the Raman- and 
IR-active modes with non-vanishing frequencies. The electron-phonon coupling 
parameter in calcium intercalated graphene trilayer is calculated to be A = 0.80. 
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TABLE III: Raman- and Infrared-active modes, with nonzero frequencies, in Ca- 
intercalated graphene trilayer. 

a;(irrep)/cm~"'^ 



Raman-active modes 


AQ{E,) 


83(^13) 


102(^3) 


249(^13) 


411(i?,) 




489(^3) 


492 (Sg) 


1112(^5) 


1131(^3) 


1133(^3) 




1326(.4ig) 


1344(^13) 


l2,m{Eg) 


1382(^1^) 


1440(^3) 


Infrared-active modes 


65(^„) 


im{Eu) 


159(^2«) 


285(A2„) 


425(^„) 




489(S„) 


4Q2{Eu) 


1111(^«) 


1132(^„) 


1133(E„) 




1326(^2«) 


1344(A2«) 


1362(^„) 


1401(^2«) 





The calculated phonon logarithmic average frequency is 25.2 meV, and the cal- 
culated superconducting transition temperature is 10.1 K. 

IV. CONCLUSIONS 

We have carried out energy band calculations, within density functional the- 
ory, of a graphene bilayer intercalated with the alkaline earth atoms A=Ca, Sr, 
or Ba, and a graphene trilayer intercalated with Ca. In all cases, we found that 
the Fermi level is crossed by an s-band derived from the intercalant states as well 
as graphitic vr-bands, in similarity to the situation in superconducting graphite 
intercalation compounds. The electron-phonon coupling parameter A is largest 
for the case A=Ca, where A = 0.60 for the intercalated bilayer, and A = 0.80 
for the intercalated trilayer; this is to be compared with A = 0.83 for the case of 
CaCg GIC On the basis of these calculations, we expect that graphene bilayers 
and trilayers, intercalated with calcium, will exhibit two-dimensional supercon- 
ductivity. Even though there are no experimental reports yet of the successful 
intercalation of calcium into a system of a few graphene layers, it is hoped that 
this work will stimulate further experimental studies in that direction 
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